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1 Introduction 

Entanglement is one of the most fundamental features of quan- 
tum mechanics. Besides its fascinating conceptual aspect it also 
plays an important role in quantum information science because 
entanglement of qubits is the essential requirement for quantum 
computing. Various systems have been considered as qubits |1]J2], 
among them the solid state ones seem to be very promising. In 
particular superconducting flux qubit has been developed in su- 
perconducting ring with Josephson Junction PH]. The junction 
playing the role of the tunneling barrier can be replaced by a su- 
perconducting quantum wire which allows for quantum phase slip 
|[5] . Recently a flux qubit based on semiconducting quantum ring 
with a controllable barrier has been proposed [6] . In this context 
the problem of entanglement of two (or more) solid state qubits is 
of great importance. It has been investigated for superconduct- 
ing flux qubits interacting via the mutual inductance, via the 
connecting loop with Josephson Junction and via the LC circuit 
[TlIHljSUlj . It was found [7j that entangled states do not decohere 
faster than the uncoupled states. This is remarkable considering 
the expectation that spatially extended entangled states could be 
very susceptible to decoherence. 
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In this paper we want to study the entanglement of distant flux 
qubits by swapping. Presented model considerations may be ap- 
plied both to superconducting or semiconducting flux qubits. 
We investigate two independently evolving subsystems each com- 
posed of a qubit exposed to a single mode of quantized electro- 
magnetic field (Fig. [T]). Contrary to the previous studies where 
the so called external approximation was used |[ID] in this paper 
we take into account the full qubit-field interaction. 

The entanglement swapping was originally proposed for photons 
[TT| and has been investigated both theoretically and experimen- 



tally P^^IITB] . Recently this idea has been used to demonstrate 
the entanglement of two single atom quantum bits each sponta- 
neously emitting a photon [H] . In our paper we use this idea to 
entangle sohd state qubits which seem to be the most scalable 
and integrable [IS] . The process of entanglement can be described 



in this case by the interaction Hamiltonian with controllable pa- 
rameters. The use of solid state qubits instead of atomic qubits 
described in |14j allows to build systems operating at microwave 



rather than optical frequencies. 

The scheme of entanglement swapping for the discussed system 
is presented in Fig. [H Each qubit Q interacts with an electromag- 




BSM 



Fig. 1. Entanglement swapping scheme 

netic field mode R leading to an entangled qubit-field (QR) state. 
This effect has been observed in a series of experiments [IB] . The 
two {QR)i {i = 1, 2) systems do not interact with each other and 
therefore the state of the whole system is a product state. If one 
then performs the Bell State Measurement (BSM) on Ri and R2, 
the partner subsystems Qi and Q2 will collapse to an entangled 
state although they have never physically interacted. To enhance 
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the qubit-field interaction the qubits can be placed into the quan- 
tum cavity. The photons can escape from the cavity e.g. through 
a less reflecting mirror P3IIT7] . To quantify the entanglement we 
calculate the negativity, we discuss the results for two different 
initial states. 

In chapters II- IV we investigate the behaviour at short time scales 
where the decoherence effects are negligible, the influence of de- 
coherence is studied in chapter V. 



2 The qubit-cavity system 



To show the idea we consider the rf-SQUID qubit ^ in the pres- 
ence of static magnetic flux 0*^^ The Hamiltonian Hq of such 
qubit can be written in a pseudo-spin notation 

= --j^^zCTz - -j^B^Gx- (1) 

We operate at T <C B^/ks in order to neglect thermal fluctua- 
tions. The diagonal term Bz in ([T]) has the form 



B, = 2y6{pL-l){Y - 0) (2) 

where Pl = 27rL (/c/0o) > = 0^^^ h is the Josephson junc- 
tion critical current, B^ is the tunneling energy between the two 
potential wells. Close to = (0o = h/2e) the ring is well de- 
scribed by the quantum superpositions of two opposite persistent 
current states. 

At flrst we describe the process of entanglement of a qubit Qi{Q2) 
with a single electromagnetic fleld mode i?i(i?2)- We model the 
electromagnetic fleld of the resonant cavity as an LC resonator 
described by Hji 

Hr = hwR {o)a + ^ . (3) 
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When the qubit is exposed to the quantized electromagnetic field 
the total fiux = 0'^^ + 0?, contains the quantum part 



h 



R 



a 



(4) 



which leads to the qubit-field coupling. After some algebra we 
obtain 



Hqr = 



hhjQ 



Oz + hjjR (^a + ^ j 



-hg (a + a^) (cr^ cos 6 — sin 6) 



(5) 



where 

uq is the qubit frequency 

hiUQ 1 



2 ^ABfr + Bl 



(6) 



the "mixing angle" |[T8] is 



e = tan-i I" 
B 



cV 



(7) 



and the coupling constant g takes the form 

g = h 



3(/3l-1) 



^ JtwrCr 



(8) 



The above considerations can be equally well performed for a 
semiconducting fiux qubit [6] with 



(9) 



where Iq is the amplitude of persistent current, 0o = h/e^ 
describes the tunneling amplitude of an electron via a potential 
barrier. 

Assuming realistic values of the parameters for superconducting 
qubit e.g. ur = 2ti ■ 50GHz, Ic = 0.5/^^4, we get g = 0.2ujr. 
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To discuss the qubit-field entanglement we assume that the co- 
herent couphng overwhelms the dissipative processes (strong cou- 
pling regime). For creation and manipulation of entangled states, 
it is thus essential that both the cavity decoherence time Tr and 
the qubit decoherence time Tq are much longer than the qubit- 
cavity interaction time T^^ ~ tt/^ ~ 10^^^ s. Recently a high qual- 
ity cavities (quality factor Q/ ~ 10^ — 10^) have been built P^lfTS] , 
They have a photon storage time Tr in the range 0.3/is — 300/is. 
The estimated decoherence times Tq of the considered qubits are 
of the order of a few fis (to be specific we assume Tq ~ Ifis [4j). 
In the next two chapters we investigate the system at t <^ Tq^Tji 
allowing the entanglement to be obtained before the relaxation 
processes set in. 

3 Entanglement swapping 

The {QR)i, i = 1,2 system is described by a state vector \il^QR{t))-, 
which at ^ = is a direct product of the qubit and the cavity 
states: 

pm)M=\^QRmn{^QRm. 

|^gi?(0))z = \(Jn)i = \a)i \n)i, (10) 

where a represents the qubit pseudo-spin states ((/-ground ,e- 
excited), \n) are the photon number eigenstates, forming the so 
called Fock basis, n = 0, 1, 2, .... 

The interaction of the qubit with the field leads, in general, to 
the entangled state 

\^QRm = e~'^^^-^\^JjQn{0)), (11) 

As the two qubit-boson subsystems do not interact with each 
other their time evolved state remains separable: 
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The time evolution of this composite is a product of two unitary 
evolutions of its constituents generated by the Hamiltonian (IS) 
where 

\'^QR{t))i = T.[anit)\gn)i + bnit)\en)i] (13) 

n 

\i^QR{t))2 = Y.[an{t)\gn)2 + 6.(^)1 en) 2] (14) 

n 

The BSM is performed on electromagnetic modes in Fock basis 
(one photon with the vacuum) [il2j and projects the formerly 
independent qubits onto an entangled state 



PQQ{t) = Trn{\Bl,){B],\p{t)), (15) 

where 

|Bi) = -J=(|01)-|10)) (16) 

is one of the Bell states of the electromagnetic field modes, the 
trace Trji is taken with respect to photonic degrees of freedom. 
After the BSM, the final qubit-qubit (QQ) state is of the form 

IV^Qq) = [ao{t)di{t) - ai{t)ao{t)]\gg) 

+ [ao{t)h{t) - ai{t)bo{t)]\ge) ^^^^ 
+ [bo{t)ai{t) - bi{t)do{t)]\eg) 
+ [bo{t)bi{t) - bi{t)bomee) 



We quantify the entanglement by the negativity |[T9] N{p) = 
max(0, — Ai), where are negative eigenvalues of the par- 
tially transposed [20] density matrix of the two qubits. For an 
entangled state, the negativity is positive reaching its maximal 
value = 0.5 for maximally entangled pure state. It vanishes 
for disentangled states. Moreover, as it is an entanglement mono- 
tone it can be used to quantify the degree of entanglement. The 
use of negativity, instead of some entropic criteria as e.g. linear 
entropy, allows for simultaneous treatment of the entanglement 
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of pure and mixed states. Let us notice that in general (e.g. be- 
yond Jaynes-Cummings approximation) the qubit-resonator sys- 
tem evolves in an infinite dimensional Hilbert space. It is known 
[2TI that in high dimensional systems the so called PPT (positive 



with respect to partial transposition) entangled states can occur. 
They cannot be detected by the Peres criterion and negativity. 
In this paper we limit our attention to the NPT entangled states 
i.e those which are negative with respect to partial transposition. 



4 Numerical results 



We present results for entanglement of both qubit-field (A^^) and 
qubit-qubit (Nqq) systems. As the calculations are numerical we 
are not limited to the weak coupling regime. In numerical cal- 
culations the Hilbert space of microwave modes is truncated at 
^max = 10. We test the validity of the truncation by controlling 
the traces of the matrices [22] being never smaller than 0.99. 



There are many parameters affecting entanglement of qubits. To 
show the idea we restrict our considerations to selected examples 
and discuss the results for two initial states. In our model calcu- 
lation we assume that both qubits are identical, the analysis can 
easily be extended beyond . In this paper we consider only the 
resonant case i.e. ur^ = loq^ = ur = 2t^ • bOGHz. The values of 
gi are in the units of lur. 

At first we assume the initial state to be 

|V^gi?(0))i |V^Qi?(0))2 = |eO)i ^ \gl)2. (18) 

In Fig. [2] we show how the qubit-field negativity depends on the 
coupling strength g and in Fig. B] its behaviour for different values 
of the mixing angle 6. Comparing these figures we see that both 6 
and g infiuence the effective qubit-field interaction strength. The 
increase of g causes the increase of the Rabi oscillation frequency 
and the entanglement arises faster than for weaker coupling. Sim- 
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Fig. 2. The qubit-field negativity for different values of g, 9 = Tr/2, initial state |eO) 
and ujQ = lor = 27r • 50GHz. 




Fig. 3. The qubit-field negativity for different values oi 0, g = 0.2, initial state |eO), 
= = 2'/r • 50GHz. 

ilarly, bringing 9 closer to 7r/2 increases the Rabi frequency. For 
^ = the QR entanglement disappears. In the following we as- 
sume 6 = 7t/2 which gives the strongest effective coupling with 
fixed g. 

In the upper panel of Fig. [4] the oscillating qubit-field negativi- 
ties A^i and N2 reflect the varying degree of entanglement as a 
function of time. The differences in these two curves arise from 
different initial states for {QR)i systems ( |eO)i \g^)2)- If we 
perform the BSM at certain time t we obtain an entanglement 
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5 10 15 20 

Fig. 4. Negativities A''i (crosses), N2 (dashed line), Nqq (solid line) and probabil- 
ities for finding the two qubits in different states after the BSM. The initial state 
|eO)i (8) 151)2, coupling strength gi = 0.2, Oi = 7r/2. 

of qubits (solid line) conditioned by the degree of entanglement 
of {QR)i. In particular if we do the BSM at the time window in 
which the {QR)i subsystems are almost maximally entangled we 
obtain the maximally entangled qubits with Nqq ~ 0.5. On the 
other hand if we perform the BSM in the time window where the 
{QR)i subsystems are weakly entangled the QQ entanglement is 
vanishingly small. 

We emphasize that the 'time' in the figures is either the physical 
time of the quantum evolution of the QR system or the time, 
called the 'BSM time', at which the BSM was performed. 

The bottom parts of Fig. [H and Fig. O show the probabilities of 
finding the qubits in |ee), \eg), \ge) and \gg) states (e.g. Peg = 
\{eg\ipQQ{t))\^) . We see that the final state belongs to the sub- 
space spanned by \eg) and \ge). This is because of the value of 
6i = 7r/2 and the chosen projection operator. For such 6 the inter- 
action term in reduces to the form g (a^ + aj that excites 
only I en) with n even and \gm) with m odd if we start from 
|eO) and l^^l) initial states respectively. Then when the BSM is 
done the only nonzero elements, in equation ([T7f). are bido and 
ai6o- The relation between the 'occupation probabilities' can be 
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directly translated into the entanglement of the state: the more 
one of the probabilities dominates the other the less entangled is 



the state and when the probabilities 
entanglement reaches its maximal value. 



and Pqe equal 0.5 the 



The decay rate of the QR system can be estimated as 




(19) 



Assuming the cavity with Qf = 10^ we find Tr ~ 0.3/iS and 
Tqr ~ 0.5/is. For the cavity with Qf = 10^ we get Tr ~ 3/is and 
Tqr ~ 1.5/is. The decoherence time of the QQ entangled state 
is accordingly Tqq ^ Tq ^ l^s. This estimation is in agreement 
with the experimental findings [7] that entangled states do not 
decohere faster than uncoupled systems. 
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Fig. 5. A^i (crosses), N2 (dashed line) and Nqq (solid line) negativities (top) and 
probabilities (bottom) for the initial state (20), 9i = 7r/2 and small detuning 
51 = 0.2, 52 = 0.202. 

For the initial state 



\i^QR{0))i \ijQR{0))2 = |eO)i (g) |e0)2 



(20) 



the situation looks different. The identity of the systems (the 
same parameters and initial states) leads to the striking results. 
Whenever we perform the BSM we almost always (with some ex- 
ceptions) obtain the maximally entangled qubit-qubit state. In 
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order to show some subtleties we take the systems shghtly de- 
tuned with gi = 0.2,^2 = 0.202 and treat g\ = §2 = 0.2 as a 
limiting case. Because the two QR systems are almost identical 
differing minutely in Rabi frequencies, they evolve to almost the 
same quantum states and even if QRs are not strongly entangled 
the BSM gives nearly the same probabilities Peg = Pge ~ 0.5. In 
consequence, we get almost maximally entangled QQ state for 
arbitrary BSM time, except for some moments (in Fig. O for 
ujRt ~ 16 ) at which the norm of the BSM output approaches 
zero and the above quantities become undefined. If the BSM were 
performed at these moments the entanglement would be unsuc- 
cessful. In the case g\ = §2 the probabilities are always the same 
and the qubits get maximally entangled for each BSM time (see 
B line in Fig. [71) with the exceptions described above. Similar 
results we have obtained for the initial state l^^l)! ® \g'^)2- 

5 Decoherence 

Design and construction of quantum devices is always limited by 
the influence of environment. Here, instead of rigorous treatment, 
developed e.g. for pure dephasing [231I21|, we apply the commonly 
used Markovian approximation [25] and model the reduced dy- 
namics of the QR system in terms of master equation generating 
complete positive dynamics [26j. Following [tlBj we assume that 
the effect of environment can be included in terms of two inde- 
pendent Lindblad terms: 



pit) = [Lh 




2^2]p(t) 



(21) 



where the 'conservative part' is given by 





whereas the 'Lindblad dissipators 



11 



Lk{-) = 4Ak{-) + {■)AiAk - 2Ak{-)Al k = l,2. (23) 



are expressed in terms of creation and annihilation operators 
'weighted' by suitable hfetimes Ai = a/\/Tji and A2 = a-/^jTQ. 
To be precise we assume Tq ^ l^s^ Tji ^ 0.3fis. As the dynam- 
ics becomes non-unitary the system evolves, in general, to the 
mixed state. The BSM applied to the density operator of the 
mixed states is well defined physical operation of projection and 
reduction which can be shown to be completely positive (see Ap- 
pendix) and thus applicable to arbitrary p. In Fig. E] we show 




Fig. 6. (color online) The QQ negativity with (A) and without (B) decoherence 
for the initial state ([T8|). The parameters are 9i = 7r/2, gi = 0.2, Tr. = 0.3/xs, 

the results of the master equation simulations of the negativity 
(the line labeled by A in Fig. E]) in comparison with the calcu- 
lations which neglect decoherence (the line labeled by B) for the 
initial state (fTSf). The periodicity with the decoherence included 
is conserved. For better visibility we present the results only in 
a short time period. We see that decoherence decreases slightly 
the amplitude of the oscillations. 

The influence of decoherence on the entanglement of the system 
starting from (1201 ) (Fig. [7j) is much more dramatic. In contrast 
to the non dissipative case (B) the result of the BSM depends 
strongly on the BSM time and the character of the entanglement 
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becomes quasi-periodic. In Fig. [HI we show the decrease of the 
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Fig. 7. (color online) The QQ negativity with (A) and without (B) decoherence for 
the initial state (j20p . The parameters are as in Fig. [6l 



amplitude of negativity as a function of time in the larger time 
scale for both initial conditions. The decrease is faster for the ini- 
tial state ( JTHf ) in comparison with that for the initial state (pUf). 
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Fig. 8. (color online) The amplitude of the qubit-qubit negativity plotted as a func- 
tion of BSM time for two different initial states. The parameters are as in Fig. 

El 
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6 Conclusions 



We investigated a mechanism for creation of entanglement of 
two qubits, each interacting with a single mode electromagnetic 
field coming from independent sources. This interaction leads to 
two independent entangled qubit-field states and that BSM per- 
formed on the electromagnetic field modes projects the qubits 
onto an entangled state. Thus we discussed transfer of quantum 
information between systems having different physical nature and 
defined in Hilbert spaces of different dimensions. 

In the first part of the paper we have dealt with the pure states 
which is justified to some extent by estimated relatively long de- 
coherence times. The discussed systems offer the advantage of 
reaching a strong coupling regime between light and matter. We 
have checked that the Jaynes-Cummings model, valid for weaker 
QR coupling [27], gives the results in agreement with our calcula- 
tions for g < O.OScJi?. Assuming reasonable values of parameters 
we found that the strong couphng regime {T^^ ^ T^^, Tq^) can 
be realized and coherent manipulations of qubits (especially with 
the quantum error correction technique) and maximally entan- 
gled qubit-qubit states are possible. 

Analyzing the dynamics of the system in the presence of deco- 
herence we found that the observation of coherent phenomena 
and in particular the generation of highly entangled states is still 
possible. 

It seems that entanglement of distant qubits by swapping can 
have some advantages over standard schemes of setting up en- 
tanglement that rely on generating entangled subsystems at a 
point and supplying them to distant The qubits emerge 

entangled despite the fact that they never interacted in the past 
and therefore they do not influence each other by disturbing the 
single qubit features. They can be at much larger distances as the 
scheme does not depend essentially on the distance between them. 
The degree of entanglement depends on the moment at which the 
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BSM was performed. Verifying experimentally that two qubits 
are unambiguously entangled is a difficult task requiring sophis- 
ticated methods such as e.g. quantum state tomography |2B]- The 
solid state qubits and their entanglement discussed in this paper 
can be scaled to a larger set of quantum bits [22]. It can be of 
interest in the study of fundamental laws of quantum mechanics 
and can be useful in quantum information processing and quan- 
tum communication. It seems that the experimental realization of 
the presented model considerations may be performed with cur- 
rently available technologies. Following [21], we hope that 'what 
is predicted by quantum formalism must occur in the laboratory'. 
Sooner or later. 
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Appendix 

We will prove that the transformation described by ( JTHI ) and ( ITBl ) 
which was defined for pure states, makes sense also for an arbi- 
trary mixed state of the two-qubit-field (QRQR) system, ie. it 
is described by a completely positive operator transforming an 
arbitrary density matrix of the full system into a density matrix 
of two-qubit (QQ) system. Although it is easy to understand on a 
purely physical basis (transformation consists of a measurement 
and a reduction to a subsystem), it is instructive to give an ex- 
plicit proof of the statement. As a bonus we will easily find an 
explicit Kraus form of the transformation in question. 
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Let 

P= E P^'^L\pkvl){,mTnl (24) 



where, cf. ( TTUI ) 



I Idkul) = \ld)i0\k)i0\ u)2 I /)2 = I i^QR)! ® I '0Qi?)2 (25) 

for fi^v G {g^ e}, A;, / G {0, 1, . . .} form a basis of pure states for 
the full system. 

For the moment let us consider only the Jaynes-Cummings ap- 
proximation where we take into account only the modes | 0) and 
I 1) of the electromagnetic field, hence all Latin indices in I HM 
and (ES) take the values 0, 1 only. In this case density matrices 
of the QRQR system act in the 16-dimensional complex space. 
Til = C^^, and as such form a subset of the 16 x 16-dimensional 
complex linear space. Analogously, density matrices of the QQ 
system, acting in the 4-dimensional complex space, li.2 = C^, 
form a subset of the 4 x 4-complex space. The transformation 
(denoted in the following by A) described by ( fTHI ) regarded on 
the whole 256-dimensional complex space transforms it into the 
16-dimensional one. Straightforward calculations give 



A{p)=:a= E <\p^){<^rl (26) 

^,Te{.9,e} 

^<^T — 2 V^<;0t1 P^OtI P?lrO ^ P^tO) j ' ) 

where | fiu) := \ fi)i <S> \ 1^)2 form a basis of pure states of the QQ 
system. In the following we will need only 



A{jfikiyl){c^mTn\^ = 



- ^OkMlm^On + SlkSoAm^On) Ip^K^^^I (28) 

To check the complete positivity of A we use the Choi- Jamiolkowski 
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isomorphism defined as [ED] 

J(A) = (A(g)Ii)(P+), (29) 

where Ii is the identity operator on the 256-dimensional space 
and P+ is a maximally entangled state on the Tii ® Hi space 

p+ = |$+)($+|, I = E I Mi^O ® I Mi^O- (30) 

According to the Choi theorem A is completely positive if 
and only if J{h) is a positive-definite operator. Applying (ESI) 
and (EHD to (EOD we get 

J(A) = |$)($|, (31) 

where 

1^) = ^El/^^) <^ - l/ilz/O)). (32) 



article Hence J{^) is a projection and as such a positive definite 
operator, consequently A is completely positive. 

The obtained results allow to write explicitly the so called Kraus 
form of A, 

m = i:AnpAi (33) 

n 

where An are dim7Y2 x dimTYi = 4 x 16 matrices. To this end [32] 
we have to perform the spectral decomposition of the positive 
definite operator J7(A) 

j(A) = Erfnix;)(x;i. (34) 

Since dn are positive, we can rescale the eigenvectors 

\Xn) ■■= \fdn\Xn)- (35) 

Now the operators An can be found in the form 

An := {h® {^+\){\Xn) (36) 

where 112 is the identity on 7-^2- The above formula should be 
properly understood. Observe that since | Xn) is an element of 
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7^2 (S) ^1, it has the form | Xn) = | <^ I ^n,i), where | G 
"^2, I ^n,i) £ "^1, whereas ($+| = J2{iikvl\ {iikvl\. Hence 

An={h(^{<^+\){\Xn)®h) 

= E {f^kul\^n,i) \(l)n,i){^kul\. (37) 

In our case J7^(A) has only one non vanishing eigenvalue corre- 
sponding to the eigenvector | ^i) = | = Y.^w \ (l>i,^iu) ® \ 
Hence | (pi^^^) = \ jiv) and | ^1,^,^) = (| - | ^lz/0))/\/2. Using 

(I32D and ( ^71 ) we obtain finally: 

^ = 4^ E (I/"^)(mOz.1| - |M(/ilz/0|) (38) 

A short calculation shows that indeed, cf. (ESj), 

A(p) = ^M^. (39) 

The calculations do not change considerably if we go beyond the 
Jaynes-Cummings approximation, by taking into account arbi- 
trary finite numbers of photons in each cavity. In fact, in this 
case, the only difference consists of extending all summations 
over the number of photons from two terms corresponding to 
and 1 to the desired numbers of cavity excitations which we 
would like to regard. The final results ( 1551 ) and (ESI) remain unal- 
tered. The situation is more subtle if we want to take into account 
the infinite number of possible photonic excitations of the cav- 
ity modes. The corresponding cavity Hilbert space becomes now 
infinite-dimensional and a straightforward generalization of the 
Choi-Jamiolkowski isomorphism does not exist - one has to re- 
sort to slightly more involved procedures to investigate directly 
the complete positivity ||33]. It is, however, not really needed in 
our case. As it is easy to check, the final result (1381). ( I39l ) is correct 
also in the infinite-dimensional setting. 
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